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Q\ t We present a model calculation of leading order interference fragmentation functions arising 

from the distribution of two hadrons produced in the same jet of a fragmenting quark in a hard 
process. Using a simple spectator model for the mechanism of the hadronization process, for 

i— -H , the first time we give a complete calculation and numerical estimates for the example of a 

proton-pion pair production with its invariant mass on the Roper resonance. 
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■ In a companion paper 111 we have discussed the general framework of two-hadron interference fragmentation 
functions, their definitions and properties. We have given the cross section of two-hadron inclusive lepton- 
nucleon scattering as an example how they can be disentangled experimentally by using azimuthal angular 
dependences of a hard process. In the present article we supplement our investigation with a spectator model 
calculation of the interference fragmentation functions for the special case of the hadron-pair in the quark jet 
1 being a pion and a proton. 
Oh, Fragmentation functions (FF), like distribution functions (DF), ultimately contain the complete information 
on the confinement of partons inside hadrons. Factorization theorems for hard processes, whenever available, 
ensure their universality, i.e. FF (and DF) are process independent; once determined in a hard process they can 
• be used to predict observables in other hard reactions. FF and DF describe the properties of partons in hadronic 
bound states at a low scale. Their intrinsic non-perturbative nature presently prevents a complete calculation 
from first principles within the framework of Quantum Chromodynamics. On the experimental side, several of 
the DF are nowadays quite well known from deep inelastic inclusive lepton-nucleon scattering, the Drell-Yan 
process or from jet production in nuclcon-nuclcon collisions. The experimental knowledge of the various FF, 
however, is very poor because it requires the ability of measuring more exclusive channels in hard processes 
(semi- inclusive cross sections, transverse momenta and/or polarizations of detected hadrons, etc.). Therefore, 
model calculations of FF are of particular importance to learn about their possible features and to estimate 
cross sections for future experiments. 

When new channels are open in the final state, several FF arise which can also be time-reversal odd [||^] (for 
brevity, "T-odd"), in the sense that the constraints due to time-reversal invariance cannot be applied. The main 
reason for that is the existence of residual Final State Interactions (FSI) between the produced hadron(s) and 
the remnants of the fragmenting quark. Thus, FF not only contain information complementary to the one given 
by DF, but in particular the T-odd ones also represent a powerful tool to explore processes taking place inside 
the jet. Moreover, specific T-odd interference FF can represent the necessary chiral odd partner to address the 
transversity distribution hi 0. Therefore, the aim of the present paper is to calculate T-odd interference FF 
within a model. 

The presence of FSI means that in the fragmentation process there are at least two leading interfering 
channels. In the previous paper |jj we have shown that the interference of two channels is not enough to 
generate "T-odd" FF. For example, in the case of one-hadron semi-inclusive processes we can model FSI 
by introducing an external potential. Despite this simplistic approach, we encounter a serious mathematical 
difficulty, because the potential in principle breaks the translational and rotational invariance of the problem. 
Simplifying assumptions about the symmetry of the potential can be introduced, but at the price of loosing any 
contribution to the "T-odd" structure of the amplitude. Alternatively, we have explicitly shown that a genuine 
difference in the Lorentz structure of the vertices describing the two interfering processes is needed to produce 
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a "T-odd" contribution jy . All this means that addressing "T-odd" FF in one-hadron semi-inclusive processes 
requires the ability of describing the FSI between the outgoing hadron and the remnant of the jet, relating the 
modifications of the hadron wave function to a realistic microscopic description of the fragments. 

Because of these arguments, a more convenient way to investigate T-odd FF is to look at the hadronization of 
a current quark where two leading hadrons are detected within the same jet, considering the remnant of the jet 
as a spectator and summing over all its possible configurations. By interference of different channels producing 
the two hadrons, FF emerge which are "T-odd", and can be both chiral even or chiral odd jy. For the case of 
the two hadrons being a pair of pions the resulting FF have been proposed to investigate the transverse spin 
dependence of fragmentation. Collins and Ladinsky || considered the interference of a scalar resonance with 
the channel of independent successive two pion production. Jaffe, Jin and Tang || proposed the interference of 
s- and p-wave production channels, where the relevant phase shifts are essentially known. We will estimate the 
FF in the case of the pair being a proton and a pion produced either through non-resonant channels or through 
the Roper (1440 MeV) resonance. 

This paper is organized as follows. In Sec. |l| a brief summary of the definitions, properties and relevant 
formulae about two-hadron FF are given together with a description of the kinematics. We restrict ourselves to 
the information necessary to keep this paper self-contained; full details can be found in Ref. jj]]. In Sec. Ill we 
discuss an extended version of the spectator model used in Ref. and here adopted for the calculation. In 

Sec. IV numerical estimates are presented and discussed for "T-odd" FF that emerge to leading order. Finally, 
a brief summary is given in Sec. [v|. 



II. FRAGMENTATION FUNCTIONS FOR PROTON-PION PAIR PRODUCTION 

As discussed in some detail in the companion paper jjj the two-hadron FF can be defined as Dirac projections 
of (partly integrated) quark-quark correlation functions. In the field theoretical description of hard processes 
those functions describe the soft parts which connect quark lines to hadrons, i.e. they are hadronic matrix 
elements of non-local operators built from quark (or gluon) fields. For a quark fragmenting into a jet which 
contains a pion-proton pair the appropriate correlation function (in a light-cone gauge) is 

A y (fe;P p ,P w ) = ^ J 0^ e*< (0|^(C)k,p,X)(X,p,7r|^.(0)|0) (1) 

where the sum runs over all the possible intermediate states containing the pair. The momenta are indicated 
in Fig. @. 
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FIG. 1. Quark-quark correlation function for the fragmentation of a quark into a proton-pion pair. 



Before proceeding with the definition of FF obtained from A, we summarize the discussion of kinematics 
as given in . For clarity we give an explicit parametrization in a reference frame where the sum of proton 
and pion momenta Ph — P p + P n has no transverse components. Two dimensionless light-like vectors n + 
and n~ satisfying n+ ■ n_ = 1 can be used to decompose a generic 4-vector a in its light-cone components 
= (a? ± a 3 ) / V2 = a.n^ and a two-dimensional transverse vector a T . Throughout the paper we use the 
notation a' 1 = [a~, a + , a T \. A possible parametrization of the momenta is 
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where P = (P p — P n )/2 is (half of) the relative momentum between the proton and the pion, and we defined 
the light-cone momentum fractions 
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With the constraint of reproducing on-shell hadrons at fixed mass (Pi = Mr, = M^), the dependence of A 
on the momenta k, P p , P n in Eq. (yj) can be reexpressed in terms of seven independent variables: the light-cone 
component of the hadron pair momentum, P^ ; the light-cone fraction of the quark momentum carried by the 
hadron pair, Zh — P t ~ /k~; the fraction of hadron pair momentum carried by each individual hadron, £; the four 
independent invariants that can be formed in this case, i.e. 
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where Eq. ([k]) can also be expressed as 



R 2 T = a (i - o mi - a - o m 2 v - i mi 



(4a) 
(4b) 
(4c) 

(5) 



In the following definitions of the FF the correlation function A will occur integrated over the (hard-scale) 
suppressed light-cone component k + at k~ = P^ / Zh and traced with a certain Dirac structure T in the form 



Ain = _L J dk+ J dk- 
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This will reduce the number of independent variables to the five %h, £, M?, a d . To make this evident we 
rewrite the integration in a covariant way using 
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We get the final expression 
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(9) 



where the dependence on the transverse quark momentum k I through ah is made explicit by means of Eqs. 
and (^|). Using the parametrization Eq. (0) it is possible to further reexpress the set of kinematical variables 
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as Zh, £, fcy an< l ■ Rt, where R T is (half of) the transverse momentum between the two hadrons in the 

considered frame. 

Two-hadron FF to leading order are defined by the following Dirac projections 

A^{z h ,£, k?,R2, k T ■ R T ) = D^zh, f , fc T 2 , i? 2 , k T ■ R T ) (10) 

A^(,a,ei 2 ,^-i) = ^M^ll' GH*h,t,k*,R*,kr ■ Br) (11) 

75 ^ (z/,,, £, k T , i? T , /c T • R T ) — ^ — H^(zfi3 j 7 ■ At) 

' ' "'' ' ■ Ht(z hl t fc 2 , i? 2 , fc T • i? T ) , (12) 



M 1 + M 2 

where we used transverse 4- vectors defined as = g^ v a„ = [0,0, a T ] (with g£ v — g^ v — n+n 1 ^ — n+nt)- In 
this manner, the FF depend on how much of the fragmenting quark momentum is carried by the hadron pair 
(zh), on the way this momentum is shared inside the pair (£), and on the "geometry" of the pair, namely on 
the relative momentum of the two hadrons (i? 2 ) and on the relative orientation between the pair plane and the 
quark jet axis (fc 2 , k T ■ R T , see also Fig. 0). 




FIG. 2. Kinematics for a fragmenting quark jet containing a pair of leading hadrons. 

In Eqs. jlc|)-(p^|), D%, G^, H^, are the interference FF that arise to leading order in 1/Q for the fragmen- 
tation of a current quark into two unpolarized hadrons inside the same jet, in the case under consideration a 
proton-pion pair. The different Dirac structures in the projections are related to different spin states of the frag- 
menting quark and lead to a nice probabilistic interpretation jlj: D\ is the probability for an unpolarized quark 
to produce a pair of unpolarized hadrons; Gi is the difference of probabilities for a longitudinally polarized 
quark with opposite chiralities to produce a pair of unpolarized hadrons; and both are differences of 
probabilities for a transversely polarized quark with opposite spins to produce a pair of unpolarized hadrons |l]] . 
Di is chiral even and "T-even" ; instead, Gj- , H^-, are (naive) "T-odd" . Gj~ is also chiral even, while H^, 
are chiral odd and, as pointed out in Ref. M, represent the natural partner in a measurement of the transversity 
distribution hi, in a a sort of "double" Collins effect [Q. 



III. SPECTATOR MODEL 



In this section we extend the formalism of the so-called diquark spectator model to calculation of the 
two-hadron interference FF, specializing it to the emission of a proton-pion pair. 

The basic idea of the spectator model is to make a specific ansatz for the spectral decomposition of the quark 
correlator by replacing the sum over the complete set of intermediate states in Eq. (|l|) with an effective spectator 
state with a definite mass Mr> and the quantum numbers of the diquark. Consequently, the correlator simplifies 
to 
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Ay [k] Pp, P-k) — 



x J (2^) 4 
0((k-P h ) + ) 



e lk < (0|^(C)|7r,p,A:)(X,p,7r|^(0)|0) 
6 {(k - P h f - M 2 D ) (0\iP i (0)\Tr,p,D)(D,p,^JO)\0) 



(2tt)3 

= Ay (fc; P p , P„ ) S(r h - u h + Ml - M 2 D ) , (13) 

where in the last line use of the definition (Q) for the kinematical invariants has been made. When inserting 
Eq. ( p^| ) into Eqs. (^0|)-(p^|), the additional 8 function allows for a drastic simplification of the Dirac projections 
(§), namely 
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(15) 



1 - Zh 1 - z h z h 

With this hypothesis, the quark decay described in Fig. [j] is specialized to the set of diagrams shown in Figs. ||, 
H and their hermitean conjugates, where the interference, necessary to produce the "T-odd" FF, takes place 
between the channel for direct production from the quark q of the proton-pion pair (p, tt) and the channel for 
the decay of the Roper resonance TZ. 
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FIG. 3. Diagonal diagrams for quark q decay into a proton p and a pion ir through a direct channel or a Roper 
resonance R. 
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FIG. 4. Interference diagrams for the same process. 

By modelling the various ingredients entering these diagrams, we can directly calculate the soft hadronic 
matrix elements of Eq. ( |l3| ) and, consequently, make quantitative predictions for various Dirac projections and 
the related FF. 

In the following, the most naive picture of the quark structure of the nucleon will be assumed, i.e. in the rest 
frame all quarks are in the l/2 + orbitals and the diquark can be in a spin singlet state (scalar diquark, indicated 
by the label S) or in a spin triplet state (axial vector diquark, indicated by A). When generally referring to 
the diquark, the label D will be kept. We will further assume that the proton-pion pair has an invariant mass 
equal to the Roper resonance one, so that the diagrams containing an intermediate Roper TZ (Figs. ||ct,^a,^c) 
will be dominant and all other channels (including diagrams §&J^c, |J> and other resonances) will be neglected. 
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A. Propagators 



Here in the following, we list the propagators needed to compute the diagrams and 
• quark with momentum k 



k 



rn 



• Roper with momentum Ph 



Ph 



\P 2 h -Ml + iM n Tn) , 



as it is quoted by the Particle Data Group (PDG) |10| for the determination of the mass and the width 
of the Roper resonance. 

(scalar/axial vector) diquark with momentum Pd = P$/Pa 
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Pd 

i ( _ Pah Pav 
P\-M\ \ M\ 

and the polarization sum for the axial vector diquark in the form y^£p ( ' A ' ) = — 9^v H ~T7T~ 

x M A 



B. Vertices 

Here in the following, we list and discuss the interaction vertices Ty needed to compute diagrams §iand|a,|c. 

The normalization coefficients have proper mass dimensions such that / d % J d 2 R T D l(Zh , £, k 2 , R 2 T , k T ■ R T ) 

be a dimensionless number to be interpreted as the probability for the hadron pair to carry a Zh fraction of the 
jet momentum and to share it in £ and 1 — £ parts. 

All vertices connecting hadron lines to parton lines will show a form factor that takes into account the 
composite structure of the hadron. A comment on the arguments of these form factors should be made at 
this point. Generally, a form factor can depend on all invariants which can be built from the momenta of the 
attached lines. In analogy with previous works on the spectator model HH we have chosen a form depending 
on one invariant only, generically denoted k 2 in the following. It is chosen as the virtuality of the external quark 
line Th whenever the corresponding vertex is attached to it, or as the virtuality of the intermediate diquark 
(quark) propagator for the An A- vertex (qDp- vertex) in Fig. |] diagram a (c). With these choices of arguments 
of the form factors, all virtual parton lines in the calculation receive an appropriate suppression preventing them 
from being far off-shell. The expressions of form factors are such that the net result is an asymptotic behaviour 
in agreement with expectations based on dimensional counting rules. 

• Roper-proton-pion (IZjnr) 

P TT 



K 



(» 



where g 2 /47r = 14.3 is the strong coupling constant of the irNN pseudoscalar interaction. Calculation of 
the Roper-proton-pion coupling fnp-K by means of two independent ways, namely the Roper decay and 
the elastic p + ir — > p + ir scattering through a virtual channel with the same mass and width of the Roper, 
produces a decay width and an empirical channel strength, respectively, that depend on fup-Tv- Comparison 
with experimental data allows fizpn to be compatible with g within the large experimental error bars. 
Hence, in the following in all vertices involving a Roper we will keep the same form and strength as in the 
ones involving a proton, also because, being the quark content the same, the asymptotic behaviour of the 
form factors should be the same. 

• quark-diquark- (Roper/proton) (qDlZ/qDp) 
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where m is the quark mass and A = 0.5 GeV is a cut-off parameter that excludes large virtualities of the 
quark. This choice of the form factor has the advantage of killing the pole of the quark propagator J|,[| . As 
remarked previously the argument of the (qDp) form factor is k 2 = t/j in diagram |(a) and k 2 = (fc-P^) 2 
in diagram |](c) . The power law a is determined consistently with the quark counting rule that drives the 
asymptotic behaviour of the FF at large Zh i.e. 



(1 " 



(1 - z h )- 3+2r+2 ^ 



(16) 



where r is the number of constituent quarks in the considered hadron, and A is the difference between the 
quark and hadron helicities. Thus, in this case a = 2. To determine the form of the vertex involving the 
axial vector diquark we took inspiration from considerations which are valid for on-shell particles, only. In 
that case, the most general form for the vertex involving the axial vector diquark is 757^ + 75(ap M + 
with a, b arbitrary parameters. The second piece can be rewritten using the difference p^ — = P% and 
the sum (p + fc) M of the momenta involved. The former is ruled out because in configuration space it is 
equivalent to taking the 4-divergence of the spin-1 diquark field, d^ipA'- the Lorentz condition d^ip^ — 
is usually taken to grant a field energy which is positive defined. The latter can be transformed, using 
a Gordon-like decomposition, into the linear combination of 757^ (to be reabsorbed into the initial first 
piece) and ^u^Pav which, contrary to the axial current, is even under G-parity transformations and, 
therefore, is forbidden [fL2|| . Hence, the general structure 757^ is assumed for this vertex. The overall 
normalizations N q s = 7.92 GeV 2 and N q A = 11.557 GeV 2 are fixed by computing the second moment of 
Di(zh) and comparing it with available data; their values, as well as the A parameter, are taken directly 
from Ref. §. 

• quark-pion-quark (qnq) 



IT 



r^qirq 



fqirq(K 2 ) = N q1T | K 2 _ A 2|a , 



where the new parameters A ff = 0.4 GeV, N q7T — 2.564 GeV have the same meaning and role as before 
and are determined in a similar manner |J. In this case, Eq. ( p^ ) gives a = 3/2. 
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• (axial) diquark-pion- (axial) diquark (AttA) 
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Only the axial diquarks are here considered because two scalar diquarks cannot couple to a pseudoscalar 
particle. The overall negative sign, together with the Dirac structure and the strength N„a, are justified 
the Appendix. The parameters A v and a are determined as before. In this case, Eq. ( |l6| ) gives a = 2. 



C. Calculation of fragmentation functions 



Using the results of the previous Sections [II A , III B we can calculate the soft hadronic matrix elements of 
Eq. ( [l3| ) corresponding to the diagrams ||aj4p^p and their hermitean conjugates; consequently, all the Dirac 
projections needed to determine the FF Di, G^v H^~, H± can be computed according to Eq. ( pT[ ) and ()lO|)-(|l^). 

In agreement with theoretical expectations the "T-even" D\ gets the dominant contribution from the 
diagonal diagram ||a, while the naive "T-odd" G^, H^, are entirely driven by the interference diagrams 
|]a,||c, thus confirming that in this case there is a close relation between the "T-odd" structure of the cross 
section to leading order and the FSI produced by the interference of the two channels for proton-pion production. 
Therefore, in the following, we will concentrate only on the "T-odd" FF. 

Before addressing the results, some remarks must be added about the flavor decomposition of FF, in particular 
about the different role played by the scalar and the axial diquarks. In fact, as demanded by Pauli principle, 
if the total spin-flavor wave function of a baryon must be symmetric, the scalar diquark must be in a flavor 
singlet state while the axial diquark in a flavor triplet state. This leads, e.g., to the known SU(4) structure of 
the proton wave function || and to the ratio 3:1 between the contributions of the scalar and axial diquarks, 
respectively, for the fragmentation of an u quark. We assume that the spin-flavor content of the Roper wave 
function is similar to the proton one. Therefore, Eqs. (0), (jlj) for the fragmentation of an u quark into a 
proton and a pion, u — > p + n, get contributions from 
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(18) 



where Aj^ a / C ^ s / A ) refers to the soft matrix element of diagram |^(a/c) involving the (scalar/axial) diquark, 
respectively. Their relative weights are given by the proper Clebsch-Gordan coefficients for the combination of 
isospins according to the above mentioned SU (4) symmetry. The latter is responsible also for the equivalence 
(u — > p) <^> (d — > n), therefore the previous equation can usefully be adopted to describe the d quark fragmen- 
tation into a neutron and a pion. On similar grounds, the fragmentation u — * n + 7r, or equivalently d — ► p + n, 
gets contribution only from the diagrams involving the axial diquark only [p| . 

By working out the matrix elements included in Eqs. ([l7|), fllq), we get the explicit expressions for the naive 
"T-odd" FF for the process u — > p + ir at leading order in the context of the diquark spectator model, i.e. 
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2 - m 2 ) 2 (77, - M 2 + 2M 2 - A/ 2 - a d ) 

x | - ( Th - m 2 ) +a d + (Af 2 - M\ + m 2 )(l - 2£) J 

^ fqnq fqSp flZpTT fqSTZ 

2 2(r h - m 2 ) 2 (77, - M 2 + 2M 2 + M 2 - 2m 2 + CTd ) 

x|(l - 20 (M 2 - AT 2 + m 2 ) (M w - M p ) + (a d + r h - m 2 ) (M n - M p ) - 2(1 - £) (77, - m 2 ) M p J 
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fqnq fqAp flZpir fqATZ 



2 Y1M\ (r h - m 2 ) 2 (77, - M 2 + 2M 2 + M\ - 2m? + a d ) 

x | - T 2 h [m(l - 20 - (2M P - M n )] -o$(m + M n ) + 2Aflm 2 [2M P (1 + - 3£M TC - m(l - 2£)] 

+ (M 2 - M 2 ) [2(M 2 - M\){m + M n )(l - 2$) - 4m 2 (£M K + m)] 

- {Mi - M\){2M p - M n + m)(l - 2£) + 2M 2 (M% - M 2 )(M p - M n )(l - 2£) 

+ 2Af 2 m 2 [£(2M P - M K ) - m(l - 2£)] 

2<r d [M p - Mr - m(l - 0] + 2M\[M n {2 - £) - M p (5 - 4£)] + 2£Af K (m 2 + Ml) 
+ (Ml - M 2 )[m(6 - 4£) + 2M W ] + m[2(Af£ + m 2 )(l - 20 - 4m£M p ] - 4£Af p Af£ 



-2a d 



M Z A [2(M P - M n ) + 0M n + m)] + m[(l - OK - m(m + £fcf w )] - M£(ZM n - M p ) 



-(M n + m)(Ml-Ml) 



(21) 



where M p , M n , Ms, Ma are the proton, pion, scalar and axial diquark masses, respectively, and M-r,, Y-r. are the 
mass and the experimental width of the Roper resonance. The explicit dependence of FF upon the invariants 
described in Eq. (||), rather than on the variables shown in Eqs. (|l7|), (p"8|), is due just to convenience of making 
the formula the most easily readable. In the following, plots will be shown of FF considered again as functions 
of Zh, £, k 2 , i? 2 , k T ■ R T , which allow for a better understanding of the underlying physics, as stressed in Sec. Q. 



IV. NUMERICAL RESULTS 

Plots of FF in Eqs. (^9|)-(pT|) are shown for the parameter values (in GeV) 

m = 0.36 , M s = 0.6 , M A = 0.8 , M p = 0.9383 , M„ = 0.1396 . (22) 



As already anticipated in Sec. Ill, the proton-pion invariant mass will be taken equal to the Roper resonance 
mass, Mh — Mn — 1.44 GeV, so that diagrams of Figs. |]&, ||c, ^6 are then negligible. The Roper resonance 
width is T n = 0.35 GeV g§. 

For the purpose of displaying our results we will choose a special kinematics where k T ■ R T = 0, namely where 
the hadron-pair plane is perpendicular to the plane containing the jet axis and the direction of the hadron pair 
Ph (see Fig. ||[ where the indicated angle a takes the value 90°). Consequently, it turns out from Eq. (Eb) 
that &d = &d(zh, £> k%, i? 2 ). Moreover, from Eq. ( [l5| ) also the invariant th depends on Zh, k 2 and from Eq. (g) 
we have i? 2 (£)- Therefore, at fixed invariant mass Mh = M-r. and k T ■ R T = the FF are actual functions of 
Zh,£,k% only. Equation (^) further constrains the variable £, because from the positivity of i? 2 (£) we have 

M^f + J^. (23) 

and, for the above fixed values of the masses, 0.43 < £ < 0.98. 

In Fig. U the interference FF tl ~* p+7r j fj^ u ^P+ 7r ; q^- u^p+n arg g k own (f rom top to bottom, respectively) 

as functions of Zh, k 2 at £ = 0.7. Similar results are obtained for different values of £ in the allowed range. The 
maximum sensitivity to the fragmentation mechanism is concentrated around the kinematical range where the 
pair takes roughly 80% of the jet energy (zh, ~ 0.8) and has a small transverse momentum with respect to the 
jet axis (fc 2 < 0.4 GeV 2 ; we recall that G±, H?, are defined as the probability difference for the unpolarized 
hadron pair to be generated from the fragmentation of a quark with parallel or transverse polarization with 
respect to its longitudinal momentum). 

It is worth noting also from the plot scales that is roughly smaller by one order of magnitude with respect 
to the other pair of FF. 
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FIG. 5. H^(zh, fcr), H 1 (zh, fcr), G~t{zh, k T ) at £ = 0.7 (from top to bottom) for the fragmentation of a quark u into 
a proton p and a pion -n. Kinematics is chosen such that the invariant mass of the pair is equal to the Roper resonance 
and k T ■ Rt = (see Fig. ^). 

By "cutting" the 3-d surfaces of Fig. || at constant values of Zh in the interesting range, e.g. for Zh > 0.6, 
one can get the trend of the dependence on transverse quark momenta in the non-perturbative low k£ region. 
In Fig. ^j, all this is shown for Gf in the same conditions and kinematics as in Fig. | (i.e., for £ = 0.7 and 
k T ■ R T = 0) and for Zh = 0.6, 0.7, 0.8, 0.9. It is evident that for increasing Zh, the fragmentation function gets 
concentrated at lower k% and have an increasingly less important tail. This result can be generalized also to 
the other FF. In other words, the more the hadron pair is leading, i.e. it takes most of the jet energy, the more 
the FF are concentrated around the jet axis with smaller transverse momentum. 

The asymptotic behaviour of the FF for very large k^ is dictated by the form factors at each vertex. It is 
instructive to look at the quark transverse momentum dependence for the contributions from each diagram 
separately. From a simple power counting in propagators and vertices (including form factors) in the diagrams, 
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k T 2 (GeV 2 ) 

FIG. 6. (?! (fex) at £ = 0.7 and zh = 0.6, 0.7, 0.8, 0.9 for the fragmentation of a quark u into a proton p and a pion 
7r. Kinematics as in previous Fig. ^ 

one can obtain the lower limit n in the l/(k£) n asymptotic decrease. Cancellations of leading power terms in 
the numerators can result in an even faster decrease of the actual results for the FF, i.e. a higher power in the 
asymptotic l/(k£) n behavior. In Table || we compare n from each diagram with the actual asymptotic behavior 
of their contributions to the FF. 





A(3aS) 


A(3aA) 


A(4aA) 


A(4cS) 


A(4cA) 


n 


3 


2 


4 


4 


3 


n{Gi) 






6 


4.5 


3.5 


n{Ht) 






5 


4.5 


3.5 


n(H?) 






5 


4.5 


3.5 



TABLE I. The lower limit ft in the l/(fcy) n asymptotic decrease as obtained from a power counting. It is compared 
with the power actually found in the FF. The power law for the contributions from each diagram are shown separately, 
where A((3/4)(a/c)(5/ J 4)) refers to the diagram (a/c) in Fig. ^|/^ involving the (scalar/axial) diquark. 



As already found in a different context |U| , the lower limit n is found to be smaller for the dominant "pure" 
channel related to the diagonal diagram ||a, i.e. it shows the slowest asymptotic fall-off in fc^ . This agrees with 
the expectation from dimensional counting rules which link the asymptotic behavior to the number of outgoing 
partons in the lowest Fock state: 3 quarks in the Roper compared to 5 (anti-) quarks in the proton-pion pair (the 
subsequent Roper decay takes place spontaneously and does not affect the parton counting in this argument). 

We have checked our results for stability under different choices for the quark mass parameter, and found a 
modest dependence. Even a drastic choice like setting the quark mass to zero results in less then 20% changes 
of the FF functions (except for the extreme end-points of the kinematical regions where the results with zero 
and non-zero m q both vanish, and quoting of a relative size does not make sense). 
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V. SUMMARY 



In this paper we have calculated the interference fragmentation functions (FF) that arise from the distribution 
of two hadrons produced in the same jet in the current fragmentation region of a hard process, e.g., in two-hadron 
inclusive lepton-nucleon scattering. 

Naive "T-odd" FF generally arise because the existence of Final State Interactions (FSI) prevents constraints 
from time-reversal invariance to be applied. This class of FF is interesting because it allows for a direct 
investigation of mechanisms for residual interactions inside jets and can contain the chiral-odd partners needed 
to isolate the presently unknown quark transversity distribution hi. 

The presence of FSI allows that in the fragmentation process there are at least two interfering competing 
channels. However, it has been shown Q that this is not enough to generate "T-odd" FF. A realistic microscopic 
description of the vertices involved in the different channels is required which naturally selects two-hadron 
emission inside the same jet as the most convenient scenario. 

To leading order, four FF arise, among which three are naive "T-odd" and are related to the fragmenting 
quark being polarized longitudinally (G^) or transversely (-fff 1 , H\)- The latter two ones are, moreover, chiral 
odd and can be identified as the partners to isolate h%. In fact, asymmetry measurements in two-hadron inclusive 
DIS on a transversely polarized nucleon target induced by an unpolarized beam can be envisaged [Q, where 
H?, enter the cross section in convolutions with h\. In particular, survives the integration over the 
transverse momentum of the fragmenting quark and can be deconvoluted from the transversity distribution. 

Calculations of such FF have been shown extending the diquark spectator model of Refs. [^-^| and specializing 
it to the case of the hadron pair being a proton and a pion produced either directly or through the Roper 
resonance. Therefore, this is the first example of an explicit, complete and detailed model calculation of "T- 
odd" FF. 

A well-known problem of this kind of approach is that it does not provide a scale dependence. In principle, it 
is a "low-scale model" , since it involves simple valence quark degrees of freedom. In Ref. a comparison with 
paramctrizations at "low hadronic scale" and available experimental data has been attempted and a reasonable 
qualitative agreement has been obtained for the known distribution functions and Di(zh). These findings 
give a good estimate for the values of the model parameters, that here are also adopted to calculate the other 
unknown fragmentation functions. Moreover, they give an indication of the level of accuracy that can be reached 
disregarding sea-quarks, gluons and evolution. 
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APPENDIX: 



In this Appendix we give a justification for the Dirac structure, sign and strength of the coupling adopted 
in Sec. [II for the (axial) diquark-pion-(axial) diquark vertex (AttA). In the following, we refer to the same 
notations and conventions of Ref. . 

The axial diquark is represented by a spin-1 field ip^ that satisfies the equation of motion subject to the 
Lorentz condition, 

P 2 ^ = m 2 ^ ; P^ = 0, (Al) 

while the field tensor is defined as = d^ip v — d^ip^. The conserved vector current is given by 

yn = i ^F'» v )*i) v - F^ip'*) , (A2) 

where F' — d^i\)' v — d u ip' and tp( ) refers to the incoming (outgoing) spin-1 field in the vertex. Analogously, 
we define the diquark axial current 

Aft = Ja.a (F'^ a - F^'*) = 2if AnA «* p p ^ - ^ P ' i/, a ) , (A3) 

whe re p' r efers to the momentum of the outgoing diquark. The Ja^a is the weak axial-diquark coupling used in 
Sec. [II B , whose sign, however, is not yet determined, and the definition of F^ v in momentum space has been 
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used. The current Aq does not fulfil the PCAC hypothesis. Analogously to the weak neutron decay, we redefine 
the axial current introducing a pole term to restore the PCAC hypothesis (or CAC, in case of vanishing pion 
mass m w ), i.e. 

A? = K + A$ = [gT - -f^)Mv , (A4) 

q z — m z 

where g M is the momentum transferred from the vertex. The pole term contains the diquark-pion-diquark (Air A) 
vertex of interest. Its explicit expression gives the Dirac structure of this vertex: 

At = - -4^A 0v = _ tSfjl^l 2if AnA e^ty/'VV- 7 - V>'*VV) 

= 4 2 g " 2 i>'* a if A .A e u0ai p v p'P ^ = 4 g " i,'* a Ti; A r . (A5) 
q z — m„ q z — ' 

In order to determine whether j a-kA is positive or negative, we need to explicitly calculate a diagram involving 
^4q , as in a hypothetical weak neutron decay n — > p + e~ + i/ e via an axial diquark depicted in Fig. |?]. Since the 
total neutron current has the V — A ^ 7^ — 7^75 = 7 M + 757^ structure, our strategy will be to calculate Aq 
from the diagram of Fig. [7] and then project out the 757^ part: the sign of the result will determine the sign of 
Ia-kA with respect to the vector current V 1 . 





FIG. 7. Weak decay of a neutron with momentum p n into a proton with momentum p p via an axial diquark. 

To calculate the diagram of Fig. |^ we use the rules described in Sec. Ill for the qAp vertex, the An A vertex 
and the axial diquark propagator. The net result is a current G M given by 

11 1 

j4, 



= / <Fk 



k 2 - m 2 {p, p -k) 2 - M\ (p n - k) 2 - M 2 
x 2if AnA fq Ap 1° tfl p ^ af>1 [{Pn - k)p g ar7 g P1 - {p P - k) p g 1<y g pa \ . (A6) 
Moreover, the Dirac-structure decomposition of a matrix M contains the 757^ projection as 

M = -Tr^M] 757 ^ .... (A7) 



Therefore, the 757^ projection of G M will be 

-±Tr[ 757 ^] = J d'k^ 



-m 2 {p p -k) 2 -M\ (p n -k) 2 -M 2 A _ (A8) 
x {~^f a, A I 2 qAp [*" (Pn +P P - 2*0" - k ■ (p n + Pp - 2k) gn} 



14 



To determine the sign of the r.h.s. of Eq. (A8), let us consider the component /j, = v = only. Noting that 
strong off-shellness of the particles is suppressed by the form factors, propagators in Eq. (|AS| ) are dominated by 
their positive energy poles. Writing the quark energy as k° = m + e with < e <C m, and taking m p = m n ~ 3m 
and Ma ~ 2m, we finally get 



|Tr [ 757 °G ] 



-8 ]AkA fqAp me 



(A9) 



Keeping in mind the V - 
and is parametrized as 



A 7 M 4- 7s7 M structure of the current, we deduce that JattA is a negative quantity 



U*a(Pa) = ~N nJ 



P\-M\ 

\n-^i\ a ' 



(A10) 



The last open problem is to determine the strength of this coupling, in other words the size of N v a- We 
can assume that the pion emission, as it would occur in a pion exchange between, e.g., two protons, is of the 
same order for all constituents of the protons, irrespective of modelling the proton as an ensemble of three 
valence quarks or of a quark and a diquark. Therefore, we can compare the amplitudes corresponding to the 
quark-pion-quark vertex, M qlTq , and to the diquark-pion-diquark vertex, MakA, assuming that they represent 
part of the diagrams describing the pion exchange be twee n the constituents (quark, diquark) of a proton and 
another particle. Using again the rules defined in Sec. Ill we can write 



\M q „ q \ 2 = Tr 



f 1 t 1 

JqTtql5 u, _ Jqirqlb- 



111 J qnq !D jl 

If — m If — m 



f 2 



(k 12 — m 2 )(fc 2 — m 2 ) 



(k-k' - m 2 ) 



\M a .a\ 2 



= Tr 



— ifAnA £ 



k a k' 



k' 2 - M 2 



f 2 

J Air J 



' \k> 2 - M 2 A ){k 2 - M 2 ) 



[(k-k') 2 -Mj_ 



k 2 - M 2 



(All) 



where k^ is the incoming (outgoing) momentum in the vertex for both quark and diquark and the couplings 
fqirqi Ia-kA are assumed pointlike for simplicity. If the virtuality of the pion is small, in the cm. of, e.g., the 
initial quark with k = (m, 0) we can assume k ■ k' ~ me, with e ~ m some energy such that e — m <C m. Similar 
arguments hold also fo r the diquark and are independent of the angular distribution of the pion. Inserting these 
relations in Eq. (All), we get 



\M a ,a\ 2 



|AV?I 2 



fq-rr 



qnq 



M A 



(A12) 



If the couplings are not pointlike and have the structure shown in Sec. IIIB, we need to search for a scale fcp at 
which they are comparable functions of k 2 , i.e. 



\fA*A{k 2 )\ = N^A 



k 2 



M\ 



|/<j7nj(fco)l fcp m 2 



A 2 I 2 



M A 



M A \k 2 -Kl\*/ 2 



(A13) 



with K = 0.4 GeV and N q7T = 2.564 GeV [§. This is true for N A7 ,a = 6 GeV and 1 < fcg < 2 GeV 2 , which is 
reasonably far from the pole A w and from the asymptotic scale where the two form factors have very different 
power laws. 
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